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A M-matrix which satisfies the Hecke algebraic relations is presented. Via the Yang-Baxterization 
approach, we obtain a unitary solution R{0, ipi, ip2) of Yang-Baxter Equation. It is shown that any 
pure two-qutrit entangled states can be generated via the universal A-matrix assisted by local 
unitary transformations. A Hamiltonian is constructed from the _R-matrix, and Berry phase of the 
Yang-Baxter system is investigated. Specifically, for ipi = <fi, the Hamiltonian can be represented 
based on three sets of SU(2) operators, and three oscillator Hamiltonians can be obtained. Under 
. - . this framework, the Berry phase can be interpreted. 
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^ ; I. INTRODUCTION 

CO . In 1984, Berry[l| predicted that a quantum system acquires a geometrical phase, in additional to a dynamical phase, 
if the environment or the Hamiltonian returns to its initial state adiabatically. Simon Q was the first to recast the 

(-H , mathematical formalism of berry phase within the elegant language of differential geometry and fibre bundles. Berry's 

O ^ theory has been generalized by extending it to non-adiabatic evolution, non-cyclic and non-unitary evolution, mixed 

^ I , \ states or non-abelian geometric phasesjp 0, S S 13- ^ ^^^t number of experiments have verified its characteristics 

^ . in various systems, including NMR [1, U, the closely related techiiique of NQR [13, [HI, ' optical systems [l3|, and 

■ so on. Recently more and more works have been attributed to it [iJI because of its possible applications to quantum 

^ . computation [isl . [l6l . \v\ . Such interest is motivated by the belief that geometric quantum gates should exhibit an 

1 intrinsic fault tolerance in the presence of some kind of external noise due to the geometric nature of the Berry phase. 
Quantum gate operations can be implemented through the geometric effects on the wave function of the systems. On 

04 ' the other hand, entanglement is a bizarre feature of quantum theory, has been recognized as an important resource 

^ \ for applications in quantum information and computation processing [ip. [20l [2ll [22il . 

. Very recently, braiding operators and Yang-Baxter e q uati on ( YBE) 23? 24 |25| h ave been introduced to the field 



o 

ON 



' of quantum information and quantum computation [26l . [27l . [2a ! [29l IsoT 31 ITiflal, [H, [Si- Ref. [1^ investigated 

o 



quantum computation by anyons based on quantum braids. Ref. [27[ has explored the role of unitary braiding 
operators in quantum computation. It is shown that the braid matrix can be identified as the universal quantum gate 
[27l . l28l . [2^ . This motivates a novel way to study quantum entanglement and the Berry phase based on the theory 
of braiding operators, as well as YBE. The first step along this direction is initiated by Zhang, Kauffman and Ge 
\ [2^. In Ref. [2^, the Bell matrix generating two-qubit entangled states has been recognized to be a unitary braid 
transformation. Later on, an approach to describe Grecnberger-Hornc-Zeilingcr (GHZ) states, or N-qubit entangled 
states based on the theory of unitary braid representations has been presented in Ref. (30| . Chen and his co-workers 
. [sil . [3^ utilized unitary braiding operators to realize entanglement swapping and generate the GHZ states, as well as 
J-j ] the linear cluster states. With the unitary 0) matrix, the authors constructed a Hamiltonian. The Berry phase 
and quantum criticality of the Yang-Baxter system have been explored accordingly. In a very recent work |33| , it is 
found that any pure two-qudit entangled state can be achieved by a universal Yang-Baxter Matrix assisted by local 
unitary transformations. However, The unitary solution ii{x) of YBE, which is presented in Ref. [ssj . only depends 
on a parameter d which is time- independent. We can't explore the evolution of the Yang-Baxter system. Motivated 
by this, in this paper we present a unitary solution of YBE, (pi, (^2), where ip\ and (^2 are time-dependent, while 
d is time-independent. Thus, we can explore the evolution of the Yang-Baxter system by constructing a Hamiltonian 
from the matrix (pi, (/?2)- This in turn allows us to investigate the Berry phase in the entanglement space. 

The paper is organized as follows: In Sec|lTl we present a 9 x 9 Yang-Baxter matrix R(0., (fi, (p2)- In the following, 
we investigate the entanglement properties of it. It is shown that the arbitrary degree of entanglement for two-qutrit 
entangled states can be generated via the unitary 7?-matrix acting on the standard basis. In fact, we can prove that 
this unitary matrix R is local equivalent to the solution in Ref.[33|. So we can say that all pure two-qutrit entangled 
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states can be generated via the universal ^-matrix assisted by local unitary transformations. In Sec lIIIi we construct 
a Hamiltonian from the unitary _R- matrix. The Berry phase of this system is investigated. Specifically, for ipi = (p2j 
the three nonzero Hamiltonian subsystems all are shown to be equivalent to oscillator systems of two fermions. We 
end with a summary. 



II. UNITARY SOLUTION OF YANG-BAXTER EQUATION AND ENTANGLEMENT 



As is known, The YBE is given by 

Ri{x)Ri+i{xy)Ri{y) = Ri+i{y)Ri{xy)R,+i{x), 



(1) 



where x and y are spectrum parameters. The notation Ri{x) = Ri^i^i{x) is used, Ri^i^i[x) implies li I2 <8) I3 ■ ■ ■ (8) 
Ri.i+i{x) <E) ■ ■ - ^ In, Ij represents the unit matrix of the j-th particle, and x = e*^ is a parameter related to the degree 
of entanglement. Let the unitary Yang-Baxter 7?-matrix for two qutrits be the form 



R,ix) = p{x)[U + Gix)Mi_ 



(2) 



where p{x) and G{x) are some functions needed to determine later on, 1^ = 1^01.;+!, and the Hermitian matrices M^'s 
(i.e., Mi = M^) satisfy the Hecke algebraic relations: MiAIi^iMi + gMi = Mij^x^i^'h+i + 9^h+\^ = aMi +/31i. 



with /3 



2a2 while a ^ 0. Substituting Eq. © into Eq. H]), one has G{x) + G(y) + aG{x)G{y) = [1 



gG{x)G{y)]G{xy). The initial condition Ri{\) = li leads to G{\) = 0, p(l) = 1. In addition, the unitary condition 
R+{x) = R~\x) = Ri{x~^) yields G{x) + G{x-^) + aG{x)G{x-^) = 0, p{x)pix-'^)[l + PG{x)G{x-'^)] ^ 1. For 
convenience, we restrict ourselves on a = 1, and /3 = g = 2. As a result, one has 

In this work, we choose basis {|00), |01), |02), |10), |12), |20), |21), |22)} as the standard basis. The 9x9 matrix 
M is realized as. 
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where qi = e^'^^ and 92 — e*'^^ , with the parameters ipi and ip2 both are real. 

The Gell-mann matrices Xu satisfy I^] = i,fxfii^Iij{X, 1^,1^ = 1, • • •, 8), where = ^A^. To the later convenience, 
we denote I\ by. I± ~ Ii i: i/2, V± = /4 =F il5,U± = /g ± j/7, Y — -^Ig- Introducing three sets of SU{3) realizations 



I. 



(1) 
± 



/^2^ c/^^ = c/^^/2^ vi'^^v,^v,^ 



l'^' = Wl + ll) + h{llY2 + Y,lll 



yd) = i(Yi + Y2: 



(4) 



{li) : < 



j(2) _ TT±^/± ttC^) 



r(2) 



2 ' 



V_ 



(2) 



^1 ^2 ' 



iHiil + ^2) + ^2) + 1IY2 + Yi/f], 



(5) 
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{^^^) : <^ if} ^ i [_ i (/3 + j3) _ i (y^ _ y^) + ifY^ + Fi/f], 



(6) 



where 4^=) = I^^^ ± i/f \ ^ = 4^=^ T i/f'\c/i''^ = ' ± ^/f \ ^t^') = -^I^Hk = 1,2,3), which all satisfy the 



commutation relation = iSij f\fj_^Iu\x, fi, v = 1, • • = 1,2,3). 

For i-th and {i + l)-th lattices, in terms of the above operators @ ([5]) (O, A/ can be expressed as follows. 



91 92 92 



<hj{2) ^ 92^(2) 

92 ^ 9i " 



92 



91^(3) + 92^(3) +^^^|3) 



92 9i 
92 91 



92 91 

We eventually arrive at the unitary Yang-Baxter matrix for two qutrits as, 



9? 



qiU 



(2) 



qiU 



(3) 



(7) 



i?,(x) = i(61, + aM,), (8) 

where a = — x and 6 = 2a: + . So the whole tensor space ® is completely decomposed into three 
subspaces. i.e. C'^ = C'^ C'^ C^. In addition, each block of ^-matrix can be represented by fundamental 
representations of SU(3) algebra. 

When one acts Ri{x) directly on the separable state |00), he yields the following family of states 



lV')yB = i(6|00) + ^|ll) + 4|22)). (9) 

J 92 92 



By means of negativity, we study these entangled states. The negativity for two qutrits is given by, 

^{p) ^ , (10) 

where denotes the trace norm of p^-^ , which denotes the partial transpose of the bipartite state p [3^. In 

fact, ^(p) corresponds to the absolute value of the sum of negative eigenvalues of p^^, and negativity vanishes for 
unentangled states [37| . By calculation, we can obtain the negativity of the state \4')yb as 

4 



J^{e) = -{sm^6+ Isine'lVl + 8cos2 6). (11) 
9 



When X = e'a , the state \ip)YB becomes the maximally entangled state of two qutrits as |?/')yb = "^('^' ^ + 

4|11) + ^|22)). In general, if one acts the unitary R(x) on the basis 100), 101), 102), llO), 111), 112), 120), 121), 122), he 
will obtain the same negativity as Eq. (fTTI) . It is easy to check that the negativity ranges from to 1 when the 
parameter 9 runs from to tt. But for 9 e [0, tt], the entanglement is not a monotonic function of 9. And when 
X = e'"3, one will generate nine complete and orthogonal maximally entangled states of two qutrits. 

It is worth mention that the entanglement doesn't depend on the parameters ipi and ip2- So one can verify that the 
parameters tpi and ip2 may be absorbed into a local operation. Actually, we can introduce a local unitary operation P 
whose form is P = diag{^, l,gi}. By means of this local transformation {PiS)P)R{9,(pi,ip2){P~^ ^P^^) = R{9), we 

can say the unitary 7?(0, (pi, </J2)-matrix ([8]) is local equivalent to the universal 7?(^)-matrix, which is the solution of 
n = 3 in Ref. [s^, where the proof of universality for a. x Yang-Baxter matrix has been presented. So the same 
as the property of ^(0)-matrix in Ref . [ssj . we can also say all pure entangled states of two 3-dimensional quantum 
systems {i.e., two qutrits) can be generated from an initial separable state via the universal R{9, ipi, (p2)-matrix ([8]) if 
one is assisted by local unitary transformations. 
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III. HAMILTONIAN AND BERRY PHASE 

A Hamiltonian of the Yang-Baxter system can be constructed from the ^(0, i^i, (/52)-matrix. As is shown in Ref. 
[3l| . the Hamiltonian is obtained through the Schrodinger evolution of the entangled states. Let the parameters 6 
be time-independent and ipi = Wit be time-dependent, where Wi = n,ja;(i=l,2;^ is a fraction in lowest terms), the 
Hamiltonian reads, 



dt 



(12) 



k=l 



where the superscript k denotes the fc-th subsystem. The k-ih. subsystem's Hamiltonian H'^^^ can be expressed as 
follows. 



8y/2hLU 



smf 



/2 /2 
+ ^iib*n2 - 2nisin6')e"2»V2y|i) _ l-{ibn2 + 2ni sin 6)6^''^^^^^ 



— {ibni+2n2 sme)e^"^W]^' + —{ib*ni - 2n2 sin6')e 



12 

V2 



12 



,-2i¥3ijj(l) 



/2 

2 (m -712 ) sin 6*4^^ + ^(ni -f 7i2)sin6iy(i'] 



(13) 



smf 



[^{ibn2 ~ ib*ni)e'^'^'-'^'hX' + ^{-ib*n2 + ibni)e-'^'^'~'<"''> I 



/2 /2 
-h^(2nisin6'- i6*ri2)e*'^=yj*^ -f ^(2ni sine* + i6n2)e-"^^y 



/2 

+ + 2^2 sin6l)e-*'^i[/|'^ + ^(-i6*ni 2^2 sine)e"^'U':'^ 



12 



12 



V2, 



V2, 



+ — {ni ~ 712) sin 6*4*' - —{ni 7^2) sin 61 Y^*)] 



(14) 



Hereafter the superscript i{i = 2, 3) denotes the second or the third subsystem. 



Based on the operators I^\x ~ 1,2,- •• ,8;fc = 1,2,3), The fc-th subsystem's Hamiltonian can be rewritten as 



follows, 



(15) 



A=l 



where C(l) = ^sVM^sine ^ _ 4V2n^ sin e - ^ 3,3). By comparing EqdTS]), with EqdTS]), one can 



(k) 

obtain as follows. 



b[^^ = ^{ni — 712) cos0sin2((pi — ip2) + ^{ni + 712) sin0cos2(iy9i — (^2), 
62^'' = —-^(711 — 712) cos cos 2((y9i — (^2) + "^("-i + n-2) sin0sin2((^i — (^2), 



B 



(1) 



72/ 



-^(7ii — 772) sin0. 



-^r' = ^7i2sin6'cos(2(p2) + ^772 cos 6* sin(2(^2) - ^77i sin6' cos(2(p2), 
^5^^ -^772 sin6'sin(2(p2) + ^772 cos6'cos(2.^2) + ^77i sin6'sin(2.^2), 
i?g^^ = ^77i sin0cos(2(^i) -f -^771 cos0sin(2v3i) — ^772 sin0 cos(2v3i), 



^77i sin0sin(2(pi) + ^77i cos0cos(2i^i) + ^772 sin6'sin(2i^i). 



. = ^(ni+n2)sin0. 



(16) 
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= ^("1 - n2)sin((^i - 


- (p2)cos6'- ^(ni 


-f 7^2) cos((^i - (^2) sin 




= ^(ni - n2)cos((pi 


- 1^2) cos 6*+ ^(ni 


+ 772) sin((/3i — (^2) sin 




~ ^(n\ — n-i) sin0, 








= —^712 COS (^2 sin6' -f 


- -^?^2 sin iy92 cos + 


-^Tii cos </J2 sin 6*, 




= —-^712 sin(y92 sin0 — 


^712 COS iy92 COS Q + 


-^771 sin (fi2 sin0, 


^6 


= — -^rii cos(^i sin6' 4 


' -^«,i sin (y9i cos + 


-^7i2 cos sin 0, 




= — -^rii sin<y9i siii0 — 


^Tll C0Siy9i COS0 + 


-^712 sinipi sin0, 


^8 


= -^(ni +712) sill 61. 







(17) 



The Hamiltonian of the fc-th subsystem, iJ(B(t)^'"'^)*^'^-', depends on the parameters B^^^ (\ ~ 1,2, - • • ,8), which are 
the components of vector B'-'^-'. Namely, B'-'^^ arc a series of time- varying parameters controUing the fc-th subsystem's 
Hamiltonian. After a periods T^*''\ Hamiltonian returns to its orig inal form, i.e. ff(B(0))W = H{^{T^^^)) . 
According to this, one can easily verify the periods of the subsystems are, T^^^ = ^ and T^*' = ^ when ^ = m{i.e. 
77,1 = 771, 772 = 1, whcrc m is integer and tjt, 7^ 0). 

To the later convenience, we denote -^77^ — 771772 + = 77, 377i + 2\/2nsm9 = a+, 3ni + ^/2ib*n = 377i — 
2V277sin6' = a_, 3ni — \/2ib*n = (3^, and 8772 + 2\/2nsm6 = S+, 8772 -I- \/2ibn = 77+, 8772 — 2\/2nsm9 = (5_, 
3772 — \/2ibn ~ 7y_ . Afi''^'s{a = +, 0, — ) arc normalization coefficients. M'^^ = 1277^ [6?7^ -I- 2\/2(ni — 2772)77 sin0 — Sti^], 
Af^^'^ = 2772, = 127i2[6n2 - 2V2(ni - 2712)77 sin 6* - 877^]; 7V|'^ = 12772[8n2 - 2V2(77i + 712)77 sin 6* -|- 8771772], 

A^o^'^ = 2772, A^i'^ = 12772(8772 + 2\/2(?7i + 772) 77 sin 6' -I- 8771772]. The eigenstates of the first subsystem are found to be 



' IS^l)) = ^i=(((77l- 2772)a++ 6772)e2*'^^|00)-f(772/3_+V27fe7777l)e2''^l|ll) + (77ia+- 772/3+)|22)), 



l4'^> = ^=(-77ie2''^^|00)+ 772e2*'^l|ll)-H (771 - 772)122)), 



l^^i^^) = ^^((("1 - 2772)a_ -h 677|)e2''¥'=|00) -t- (772/?+ - ^/2i677nl)e2*'^l |11) + (ma^ - 772/3_)|22)), 



(18) 



with the corresponding eigenvalues E^^ = ^^hnuj siii6, Eq'^^ ~ and i^i^' = —^^hnujsniQ. For the second and the 
third subsystems, the eigenstates are found to be. 



E_ 



(2) 



E^ 



(2) 



E, 



(2) 



E_ 



(3) 



E^ 



(3) 



E, 



(3) 



^((77ia_ 772(S_)e*'^^|01) + (77ia_ - 772/3^)|12) -f (772(5- - 77 i77+)e-*('^i-'^^)|20)), 



^((-771 + 772)e*'^^|01) + 77l|12) - 772e-*('^l "'^^ > 1 20) ) , 



^((77ia+ + 772<5+)e''^^|01) + (77ia+ - 772/3;)|12) + (772(5+ - 7717/- )e-'('^i -'^^) 1 20) ) , 



^((772(5- - 77l?7+)e-*('^l-'^=)|02) + (77ia_ 772(5_)e*'^^ |10) + (77ia_ - 772/31)|21)), 



1 (_„2e-*('^i-'^=)|02) + (-771 + 772)e"^=|10) +77i|21)), 

^((772(5+ - 77i77_)e-^('^i-'^-)|02) -f (77ia+ + 77 2(5+ )e^'^- 1 10) + (77ia+ - 772/3;)|21)), 



(19) 



(20) 
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with the corresponding eigenvalues E^f' = ^^hnuj sinO, Eq^ = and E^^ = —^^hnu siiiO {i = 2,3). 

According to the definition of Berry phase [H, when the parameter B'-'^'' is slowly changed around a circuit, then 
at the end of circuit, the cigenstates \E^''){a = +,0,— )evolves adiabatically from to T^'^^ the Berry phases 
accumulated by the states \Ea^) are, 



(21) 



By substituting ([T5|) ([TO)) ([^0]) into (pij) . one can obtain the Berry phases for these eigenstates. 



Subsystem 1 



(1) 



ifi+K.2^r(i) 



2n2 



^(2) ^^-^ 



Subsystem 2 or 3 : < 



7o 



-i2(ni-n2)(ni-2n2) 
2n2 



(22) 



where 



A'4 



-lOn^ + 13njn2 + lining - %2n\n\ + Mn^n 
-6V2 sin 9nn2{nl - Qn\n2 + 12?! 



!cos(26')(ni - 2712)71", 



l"-2 



871 



2^7 



2n'^{Qn\{ni -~ 712) - 8sin^6'(7i? + nl)) - 'dn2{n\ - n\n2 + 57ii7i2 - 37ii7i^ + 2713), 
GV^ sin 07171,2 (77^ + &n\n2 — 2inin\ + 4772). 



The above Berry phases are the general result for the whole system. Despite the complexity in the expressions for 
the Berry phase in its full generality, we can use the general formulae to discuss some special cases for convenience 
to understand. 

Example 1: For ipi ~ ip2, i.e. 7ii = 712 = 1. In this case, T^^^ ~ ^ and T'*' ~ When one substitutes 
the conditions into Eq. ([22]) . he then gets the following explicit expressions of the Berry phase (all phases are defined 
modulo 27r throughout this paper), 



(fe) 
7+ 

(k) 
To 



-7 



0, 



(fe) 



-7r(l 



2\/2 , 



(23) 



where k = 1,2,3. In fact, for cpi = (p2, the Hamiltonian can be represented based on three sets of SU{2) operators. 
Under this framework, the Berry phase can be interpreted. In the following, we have a thorough discussing for the 
case. 

Introducing three sets of SU(2) realizations in terms of three sets of operators (O ([5]), 





^(id^-^ 








+ [/i'^) 


oik) 

'-'3 — 


|y(fe) _|_ 





(24) 



They satisfy the algebra relations of SU{2) group: [s'f\ S^l^] = 2S,jS^^'', [S^'',S^^'] = ±,5„S'^'^ {S^^>)^ = 0{i,j 



jik)y 
:(0 



(fc) oik) 



1, 2, 3), with S'j_ = S{ ±75*2 (fc = 1,2,3). By the way, their second-order Casimir operators are J'^'^^ = ^(S'_j_ 5 



5.(fc)^Cfe)^ + (4''^)2. One can verify that the eigenvalues of j'-''^ are ^(i + 1) = | and 0(0 + 1) = which correspond 



to spin-i system and spin-0 system. In terms of the operators (|24p . the Hamiltonian for subsystems Eq. p3p and Eq. 
(fT4l) can be rewritten as follows: 



ij(i) = C{l)i~lib*e''^S^^'^ + Itbe-''^S['^ + ^smeS^'^), 



(25) 
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= C{2){-lzbe''^S^^^ + Ub*e-'^S^y + ^sin^^f ), (26) 

D D O 

ij(3) = C{3){~ltbe"^S^^^ + l:ib*e~''^S^^^ + ^sin^^^^). (27) 

D D o 

C(l), C(2) and C(3) are defined in Eq. p^ . So we can say the whole system is equivalent to three spin-i subsystems 
and three spin-0 subsystems. Actually, we can introduce a 9 x 9 orthogonal matrix P which is time-independent (see 
Appendix A). By means of P, the whole system's Hamiltonian H and Casimir operators J'^''^ are transformed into 

block-diagonal matrices. Namely, H = PHP'^ and j'-''^ = Pj'^^'^P'^ are block-diagonal matrices, where P^ denotes 
the transpose of matrix P. 

For the subsystem 1, from Eq. (|A2p we can obtain its Hamiltonian if'^' ~ H^l'^ ® h'^\ For hI^\ the eigenvalue of 

2 

Casimir operator J''-^' is 0, and the Berry Phase is 0. So we can say the subsystem Hamiltonian H^^^ is equivalent to 
a spin-0 subsystem. For H^^\, we introduce the transformation, cos a = 2^ sin 6* and cos/3 = - sin e cos 2y+3 cose sin 2y 

^ h '3 I- ^9_8sin2e 

where a £ (arccos arccos — ^y^) and /? £ [0, 27r]. a is time- independent, and (3 is time-dependent. Then the Berry 
phase Eg ([23)1 can be recast as, 

7i'^=T7r(l-cosa)=T^ (28) 

where Vl{C) — 27r(l — cos a) is the familiar solid angle enclosed by the loop on the Bloch sphere, and the parameter a 
comes from d which comes from the Yang-Baxterization of the Hermitian matrix M . So the Berry phase depends on 
the spectral parameter. By means of P , the eigenstates can be recast as follows(we neglected the global phase 

factor), 

|£;W) = _e'/3sin||00)-Hcos||01), (29) 
l^fi)) = cos||00)+e-*^sin||01). (30) 

The Hamiltonian i/*^^' in Eq (P5)) can be recast based on the operators Ea (|A3p to the form 

Ef' = -2nwcosa(2cosa5^^^ +sinae''^S'';^' +sinae~*'^S'i^') 

2 

= -2nuj cos aH^^\ (31) 
where Hq^^ is of the following form, 

i?^^' = 2cosaS';^^' -Hsinae*''^^^' -fsinae-*^5i^^ 

ViOQ find OQTYio T-\Tn7-C)i rio 1 -mooniTirr o o -J-ViQ-i" rrTT7-i::^n in f^oll T'Vio^- ic TT^^^ 



Thus, the Hamiltonian Hi has the same physical meaning as that given in [33 • That is, H\ ' is an oscillator 

Hamiltonian formed by two fermions with frequency 2ti;cosa where a S (arccos arccos — 2^). For the Yang- 
Baxter subsystem, a = is a critical point. However, because of a 7^ 0, the subsystem Hamiltonian Eq (|3ip can't 
reduce to the standard oscillator. In other words, sinae*'' plays a role of the "energy gap" and the wave function 
of the subsystem takes the form of spin-coherent state (38j. The quantum criticality can't occur in the Hamiltonian 

subsystem and the Berry phases cant't vanish since a 7^ 0. 

2 

Via the same method, the Berry phases for subsystems 2 and 3 may be obtained, 7^!'' = T''i'(l — cos a) = ^ ^^^^ 
and 7^'^ = 0(i = 2, 3). The Hamiltonian ij(') in Eq(l26|) and EqdST]) can be recast to the form 

nf ^ ^ i/f ^ = -hij cos aH^''^ , (32) 

2 2 

where H^n^ = 2cosaS'l'^ + sin ae^'^S'^ + sin ae'^^S'i'^ (i 2, 3). Thus, the Hamiltonian h[^^ and h[^^ both have the 

2 2 

same physical meaning as that given in [11. That is, ifP' and h[^'^ both are oscillator Hamiltonians formed by two 
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(2) (3) 

fcrmions with frequency wcosa. Thus, in the two Yang-Baxter subsystem H\ ' and H\ , the quantum criticahty 

2 2 

can't occur and the Berry phases cant't vanish since a 7^ 0. 

Example 2: For ipi ~ — <P2, i-e- = —n-i = —1. In this case T^^^ = ^ and T^*) = One can get the following 
explicit expressions of the Berry phase, 



f (1) _ 76 sing 2^ 



Subsystem 1 : < 



7^=0 



Subsystem 2 or 3 : < 



7+ — 3 
7« = 



(33) 



In this case, one can find that 7^^' = — 71'', 7o"'"'' = — 7o*'' and 7i"'"'' = ^1+^ ■ 

Example 3: For 951 = 2ip2, i.e. ni = 2n2 = 2. In this case T^^'> ^ ^ and T^') = 221. We then gets the following 
explicit expressions of the Berry phase. 



Subsystem 1 : < 



(1) _ Vesing p 
'+ 3 

7^^^=0 



^(1) ^ _ v^27r 



Subsystem 2 or 3 : < 



(j) _ VesinSr, 

7+ - 3 
7^'^ - 



(34) 



„,(') _ _ V6sin_£ 
3 



2n 



In this case, one can find that 7^^'' = —7!*'', 79"'"'' = ~7o''' and 7!"'"'' = ^7+''- 

Example 4: For (fi = ~2ip2, i.e. 711 = -2^2 = -2. In this case T^^) = ^ and T^') = ^. We then get the following 
explicit expressions of the Berry phase. 



Subsystem 1 : < 



7i^^ 



(I + 4^ sin 0)2^ 



2tt 
7 



(4 _VT±sme)2T: 



Subsystem 2 or 3 : < 



(0 _ 27r 
" 7 



76 
{i) 

7- 



(35) 



;| + ^sin0)27r 



In this case, one can find that 7^^'' = —7!*'', ^q^^ ~ ~7o''' ^^"^ 7-"* 



-7+ • 



IV. SUMMARY 



In this paper, we have presented a 9 x 9 Af-matrix which satisfies the Hecke algebraic relations and derived a unitary 
(/3i, (/32)-niatrix via Yang-Baxterization of the A/-niatrix. In the following, we can say that all pure two-qutrit 
entangled states can be generated via the universal _R-matrix assisted by local unitary transformations. Specifically, 
the arbitrary degree of entanglement for two-qutrit entangled states can be generated via the unitary R matrix acting 
on the standard basis. Then the evolution of the Yang-Baxter system is explored by constructing a Hamiltonian from 
the unitary ^-matrix. In addition, the Berry phase of the system is investigated, and general expressions of Berry 
phase are figured out. Finally, we use the general expressions to discuss some special cases. For Lp\ = (/'2, based on 
three sets of SU(2) operators, the Hamiltonian has been represented and the three nonzero Hamiltonian subsystems 
all are shown to be equivalent to oscillator systems of two fermions. Under this framework, the Berry phase can be 
interpreted. 



This work was supported in part by NSF of China (Grant No. 10875026). 
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APPENDIX A: BLOCK-DIAGONALIZE H AND J 



(k) 



The 9x9 orthogonal matrix P which is time-independent reads, 



P 



/ 75 



10 





10 

^ 

V ^ 



1 

\/2 



1 

















1 

t 















n 
u 


n 
u 


1 
i 




















1 




























1 













(Al) 



The orthogonal matrix P satisfies the relation PP^ — P^P — Igxg, where P^ denotes the transpose of matrix P. 
By means of P, the Hamiltonian H for (pi — ip2 can be recast as follows, 



H = PHP' 



2 2 2 



(A2) 



k=l 



where H^''^ = ijf ' ® i/f 's are 2 x 2 matrices, and H^'^^'s are 1 x 1 matrices with H^"' = (0). The (3 x 3)- 

2 2 

dimensional interactional Hamiltonian system is decomposed into six subsystems. Three sets of SU{2) realizations 
can be recast as, 



= |00)(01|, 5« = |01)(00|, 5^ = ;^(|00>(00| - |01)(01|), 



(1) 



1 



(A3) 



5(.')-|ll>(20|, ^i'' = |20)(ll|, = i(|ll)(ll|-|20)(20|). 



(A4) 



5(3) _ |0l\/00l C^^) _ lr,r,v /r,-, I o(3) 1 



5r = |21)(22|, = |22)(21|, 5^ - 2(|21)(21| - |22)(22|), (A5) 

and the seconde-order Casimir operators are j'-^^ = |(|00)(00| + |01)(01|), j'^) ^ |(|ii)(il| + |20)(20|), J^^) ^ 
|(|21)(21| + |22)(22|). 
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